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ABSTRACT 

Over  the  past  decade,  advances  in  the  areas  of  propagation  and  rough  surfac'  staticnng  have  made 
attractive  the  concept  of  generating  hybrid  models  't  investigate  the  impact  o?  realtsnc  two 
dimensional  sea  surface  roughness  on  long  range  acoustic  propagation.  The  development  of  such 
models  requires  that  efficient  and  accurate  coupling  schemes  be  created.  In  'hi.v  paper,  one  aspect  of 
this  coupling  problem  is  discussed;  generating  a  suitable  representation  of  the  mcidem  field  for 
input  into  the  scattenng  model. 

KEYV,'ORDS 

Hybrid  models;  plane-wave  spectrum;  scattering;  propagation 

INTRODUCTION 

In  the  past,  one  of  the  standard  simplifications  applied  to  the  problem  of  modeling  the  transmission 
of  sound  through  an  ocean  waveguide  has  been  to  treat  the  air/sea  interface  as  a  pcrfcciiy  fiat 
pressure-release  (Neumann  B.C.)  surface.  In  this  approximation,  the  redistribution  of  energy  that 
results  from  the  acoustic  interaction  with  the  sometimes  quite  rough  sea  surface  is  ignored  On  the 
other  hand,  researchers  interested  in  the  underwater  acoustic  scattering  from  a  rough  sea  surface 
have  generally  simplified  their  world  by  neglecting  propagation  effects  that  result  from  a  depth  and 
range -dependent  index  of  refraction.  In  the  last  decade,  there  have  been  significant  advances  in  the 
modeling  capability  for  these  two  aspects  of  the  underwater  acoustics  problem  and.  just  as 
imponant,  there  has  been  a  tremendous  increase  in  computing  power  available  to  researchers  and 
modelers.  In  fact,  much  of  the  work  that  needs  to  be  done  now  centers  around  the  creation  of  hybrid 
models  that  can  take  advantage  of  advances  in  both  of  these  modeling  areas  Such  hybnd  models 
could,  for  example,  be  used  to  explore  the  effect  of  different  sea-state  conditions  on  long-range 
propagation.  A  significant  difference  between  the  technique  pursued  in  this  paper  and  some 
previous  approaches  (Buckcr,  1970;  Head  et  al..  1989)  is  that  this  approaches  allows  for  the 
redistribution  of  the  energy  due  to  the  scattering  where  it  occurs  in  space  rather  than  treating  the 
scattering  as  a  simple  loss  mechanism  averaged  over  range. 

Hybrid  models  of  the  type  alluded  to  make  use  of  existing  propagation  and  scattering  models  and 
are  (in  the  simplest  sense)  based  on  coupling  schemes  that  perform  two  different  functions  The 
first  function  is  to  determine  from  the  propagation  model  that  pan  of  the  acoustic  field  that  interacts 
with  the  bounda^  of  interest  and  to  express  it  in  a  form  that  the  scattering  model  can  accept.  The 
second  function  is  to  coherently  couple  the  resulting  scattered  field  to  the  propagating  field.  VVTiile 
both  aspects  of  the  problem  are  equally  important  and  must  be  addressed,  in  this  paper,  the  focus  is 
on  a  numerically  efficient  technique  aimed  at  accomplishing  the  former.  In  particular,  the  technique 
described  is  appropriate  for  any  propagation  model  that  can  provide  the  complex  acoustic  pressure 
field  in  depth  specified  ranges  from  the  source  and  any  deterministic  scattering  model  that  can 
take  as  input  a  description  of  the  incident  field  in  terms  of  plane  waves.  This  technique,  which 
essentially  estimates  the  arrival  structure  or  plane  wave  spectrum  (PWS)  incident  on  a  particular 
patch  of  an  assumed  flat  ocean  surface,  provides  precisely  the  information  generally  required  as 
input  to  typical  scattering  models.  The  organization  of  this  paper  is  as  follows;  first,  the  technique 
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used  to  extract  the  PWS  from  the  propagauun  model  ss  expiauied  Next,  a  tenciinuirk  ^ 

the  total  surface  reflected  field  is  developed  ftcm  lu-rmal  rtu-vk  meofs  ruiaih  a  tiuuie'u  aJ  r>,.c:  . 

IS  ci'nsidcred  the  results  of  uhich  ailow  for  the  vcn’i,.  .men  e;  ts<.o  dotefe'H  asf'ects  i-r  U'.c  i  '»%  s 
extraction  technique;  the  angular  distribution  and  the  total  energs  ictlccied  troni  tjic  sea  si.'!t.3.c 
versus  range. 

PWS  EXTRACTJON  TECHNIQUE 

ffir  most  deienninisiic  and  stochastic  scaaenng  nuxSeiv  she  necesvar;.  input  mtoniratson  tv'  , 
consists  of  some  description  of  the  pressure  held  incufent  i.m  the  s.at'enng  patch  as  ueii  a>  s.-n  r 
determintstic  or  statistical  desenption  of  the  surface  foughness  Restrutir.g  the  d;s».uss;o!',  tc  e 
problem  of  scattering  from  a  deterministically  known  surface,  it  is  clear  that  what  the  coupling 
technique  must  provide  to  the  scatienng  mode!  is  the  magnitude  and  phase  for  each  of  the 
multipaths  incident  on  the  panicular  patch  of  interest  A  panicularly  convient  description  of  tf,;s 
arrival  structure  is  in  terms  of  its  plane  wave  spectrum  (PWS)  Here  the  PWS  is  dehned  to  cornsvi 
of  the  magnitudes  and  relative  phases  of  the  plane  waves  mctdcni  cm  a  panicular  patch  of  the  (xean 

The  approach  employed  here  to  estimate  the  PWS  from  the  solution  ot  the  one-way  wave  equation 
is  based  on  the  observation  that  the  difference  in  the  venicaJ  energy  spectral  denssises  between  !wi> 
nearby  ranges,  ri  and  rj 

has  a  clear  interpretation  in  terms  changes  in  the  vcntcal  energy  distnbution  of  the  forward 
propagated  field  or.  equivalently,  in  the  redirection  of  energy  flow  between  r;  and  rx  Here  Ptr.  st  is 
the  Discrete  Fourier  Transform  (DFT)  in  z  (As  =  i/(NAz) )  of  the  pressure  field  pir,  z).  In  a  typical 
deep  water  ocean  environment,  D(ri.  T2.  s)  will  be  dominated  by  refraction  from  within  the  water 
column  and  ocean  bottom,  reflections  and  diffractions  from  bathymetric  features,  and  reflectior. 
from  the  assumed  flat  sea  surface.  Of  course,  D(ri.  rj,  s)  will  also  include  the  effeci-s  of  ancnuaticn 
but  this  can  be  assumed  to  be  negligible  provided  ra  is  not  too  much  larger  than  rj  Of  the 
mechanisms  noted,  it  is  clear  that  plane  wave  reflections  from  the  sea  surface  has  a  distinct 
signature  in  the  transform  space  that  may  be  exploited  for  this  panicular  problem. 

Consider  a  plane  wave  having  traveling  upward  (vertical  wavenumber  kz  -  2ks)  in  a  homogeneous 
waveguide  and  reflecting  from  a  finite,  mirror-like  sea  surface  centered  between  ri  and  rj  The 
energy  redistribution  due  to  the  reflcaion  from  the  surface  (ignoring  edge  effects)  is  characterized 
by  a  net  loss  of  upward  going  energy  in  the  vcnical  energy  spcctnim  at  wavenumber  k^  and  a  equal 
gain  of  energy  at  wavenumber  -kz.  In  a  more  realistic  underwater  environment  where  the  index  of 
refraction  varies  with  depth  only,  energy  that  refracts  between  two  close  ranges  has  the 
characteristic  that  the  flow  of  energy  is  from  kz  (upward)  to  kz  -  Akz  for  downward  refraction  or  to 
kz  +  Akz  for  upward  refraction.  Thus,  it  is  clear  that  except  for  the  case  where  energy  is  incident  at 
shallow  grazing  angles  (kz  small,  kz  =  Akz),  the  energy  redistribution  associated  with  a  reflection 
from  a  flat  sea  surface  is  easily  distinguished  from  energy  travebng  in  the  upward  direction  that  is 
refracted  downward  near  the  surface  or  elsewhere  in  the  waveguide  Funhcrmorc.  the  ambiguity 
that  occurs  at  shallow  grazing  angles  can  be  minimized  simply  by  forcing  ri  and  rj  closer  together" 
As  a  practical  point,  this  ambiguity  need  not  be  resolved  entirely  since,  except  for  the  roughest  of 
sea  surfaces,  if  a  reflection  is  detected  when  none  occurs  the  effect  of  scanering  from  such 
extremely  low  grazing  angles  is  typically  insignificant. 

Finally,  it  is  noted  that  the  diffractions  and  reflections  that  result  from  bathymetric  features  can  be 
excluded  from  the  energy  spectral  density  estimates  simply  by  employing  a  suitable  depth 
dependent  filter.  A  filter  of  this  son  also  ensures  that  the  generally  upward  refraction  taking  place 
in  the  lower  half  of  the  waveguide  and  the  generally  downward  refraction  of  energy  in  the  upper 
half  of  the  waveguide  do  not  act  together  in  such  a  way  as  to  provide  the  distinct  signature 
associated  with  a  reflection  from  the  sea  surface. 

The  above  discussion  suggests  die  following  algorithm  for  extracting  the  PWS  reflected  from  the 
sea  surface  between  two  ranges  ri  and  tz: 

1)  Determine  difference  in  the  energy  spectrum  as  in  Eq,  1  after  applying  a  depth  dependent 
filter  that  cuts  on  several  acoustic  wavelen^s  below  the  sea  surface. 
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2)  For  each  positive  value  or  the  venual  wavenumt'er  luehncJ  ncie  u.  Oc  upivoiOi  icvi  it 
D<k/)  <  0,  and  D(-k/)  >(1.  uidy  these  wavenumbers quaiils  as  potcntiai  s*Jf1ai,e  terlewtioriS 

3)  Since  some  energy  imbaJance  is  bound  to  occur  due  irnpan  to  tfie  iruncatuMi  or  ihe  anass 
estimate  the  retlected  energy  to  be  the  minimum  of  Dtk/i  and  Dt  -k/  i  Nonriaii/e  trie  awociated 
ampiuude  by  dividing  by  r;  -  r.  to  get  the  intensity  of  the  plane  wave 

4)  The  grav'ing  xngie  asscKiated  w  uh  this  plane  vs  jve  is  y  =  assniti.  ’k  /  •  j  2*r.*f  L  w  here  l  o 
the  sound  speed  at  the  surface  and  f  is  the  acoustic  frequency.  Ihe  relative  phase  of  the  piane  wave 
is  estimated  to  be  the  phase  (in  that  k^  bin)  given  by  the  tianslotm  of  the  venicai  pressures  at  ri 

In  the  following  section,  a  benchmark  solution  for  the  total  reflected  field  ai  the  sea  sur!av:e  is 
developed. 

NORMAL  MODE  THEORY  FOR  THE  SURFACE  REFLECrED  FIELD 

The  acoustic  pressure  in  a  range  independent  ocean  waveguide  can  be  expressed  in  a  normal  mcvde 
representation.  Using  the  separation  of  variables  technique  to  solve  the  wave  equation,  we  obtain  a 
depth  and  range  dependent  pan.  The  normal  mode  portion  Zn<z)  (venicai  eigenfunction)  is  the 
solution  of  the  following  depth  equation 


d:[  di  j  le-jr)  j 


(2) 


where  p(z)  is  the  density  of  the  medium,  o)  is  the  angular  frequency  of  the  cyiindncally  spreading 
wave,  c(z)  is  the  compressionai  sound  speed  of  the  medium  and  Kn  is  the  honzonta)  wavenumber  or 
eigetivalue.  V'hen  assuming  that  the  waveguide  consists  of  many  isovelochy  sound  speed  layers, 
the  eigenfunciion  is  expressed  as 


Z,(z)  =  A,sin(r,r)+  a,  cos(  y,:) 


(3) 


where  yn  is  the  vertical  wavenumber /eigenvalue.  The  solution  to  the  range  equation: 


^'Rfr)  .  1  dR{r)  . 


KlRir)  =  Q 


(4) 


is  the  cylindrical  Bessel  function  of  the  first  kind  Jo(Knr)  which,  considering  only  the  outgoing 
wave  and  large  arguments,  takes  the  following  form 


R(r) 


(5) 


Since  the  depth  eigenfunctions  are  orthogonal,  (and  are  required  to  be  onhonorma!)  they  can  be 
used  to  expand  a  particular  dependence  of  pressure  on  z  as  a  sum  of  the  eigenfunctions.  For  this 
case  the  function  is  a  point  source  at  r=0,  z=zs: 


S{z-z,)=  ia.Z„(r) 


(6) 


where  Qn  are  constants  to  be  determined.  Multiplying  both  sides  by  Zm  (z)  and  integrating  yields: 


Po^I2„iz)S{^-z,)^^2  =  J.Q^a\~Z„{^)Z^{z)d^ 


O) 


v'hich  reduces  to 
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where 


finally,  we  h:ue 


Po'^-i-,) 


!0. 


The  sound  pressure  for  ihe  point  source  becomes 


P{r,z)=:  C  5.- 


where  C  is  a  complex  numerical  constant.  In  the  top  layer  of  the  ocean,  ihe  boundary  condstionv 
require  that  the  coefficients  of  the  cosine  term  in  Eq,  3  vanish.  The  remaining  sme  term  which  can 
then  be  expressed  in  terms  of  exponentials  allows  the  sum  over  normal  modes  to  be  rewritten  as  a 
sum  of  up  and  down  going  plane  waves.Thc  surface  reflected  pressure  field  is  determined  by 
throwing  away  the  downward  going  plane  waves  and  evaluating  the  remainder  at  the  surface  t?  = 
0). 


2  m 


(12) 


NUMERICAL  EXAMPLE/BENCHMARKING 

The  panicular  PE  model  employed  for  this  numerical  example  was  developed  at  NOARL  (Collins, 
1988).  Several  journal  articles  have  been  written  about  this  model  and  no  attempt  will  be  made  here 
to  review  this  large  body  of  work.  The  PE  approach  (Tappen,  1977)  is  widely  regarded  as  being 
accurate  and  efficient  for  long-range  propagation  at  low  to  moderate  frequencies.  Again,  for  the 
purposes  of  this  paper,  it  is  sufficient  to  point  out  that  the  essential  attribute  of  the  PE  model,  in 
terms  of  the  implementation  of  the  PWS  extraction  technique,  is  its  ability  to  produce  the  complex 
pressure  field  as  a  funaion  of  depth  at  specified  ranges. 

In  this  section,  a  numencaJ  experiment  is  described  and  two  independent  theories  arc  invoked  to 
benchmark  or  test  the  PWS  extraction  technique.  The  range-independent  ocean  environment 
considered  in  this  example  has  a  depth-dependent  sound  speed  profile  (the  Munk  profile)  and  a  cw 
source  operating  at  250  Hz  located  at  a  depth  of  212  m.  The  waveguide  is  30Cfe  m  deep  below 
which  is  15  m  of  sediment  followed  by  an  isovelocity  half  space  which  includes  the  standard 
attenuation  ramp.  A  contour  plot  of  transmission  loss  as  calculated  by  the  PE  model  for  this 
environment  is  shown  in  Figure  1.  Note  that  in  addition  to  the  convergence  zone  which  starts  at 
approximately  40  km,  some  strong  high  angle  energy  from  bonom  bounce  paths  appear  to  reach  the 
sea  surface  beginning  at  about  10  km  and  ending  at  about  40  km.  A  ray  trace  (not  shown)  indicates 
that  these  bottom  bounce  paths  and  other  multiple  bounce  paths  contribute  to  the  field  at  the  surface 
at  all  ranges. 
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Fig,  1  C''ntour  plot  of  transmission  loss  at  93  dB  for  benchmark  problem 

As  mentioned  earlier,  there  are  two  aspects  of  the  PE -extracted  surface-reflected  PWS  that  are 
conveniently  tested.  First,  the  total  reflected  intensity  from  any  particular  patch  of  the  ocean  surface 
can  be  calculated  from  normal  mode  theory  just  described.  In  Figure  2.  a  comparison  of  the  total 
reflected  intensity  calculated  by  summing  over  the  surftce-rcflecied  PWS  and  from  the  normal 
mode  theory  for  the  reflected  intensity  (derived  from  Eq.  12)  are  compared.  For  this  comparison, 
the  vertical  filter  cut  on  at  a  depth  of  50  m,  the  range  separation  between  the  two  arrays  (r;  -  rj) 
was  2  PE  range  steps  (10  m  in  this  case),  and  a  1024  FFT  was  performed.  It  is  clear  that  both  in 
level  and  in  character  the  two  predictions,  panicularly  from  20  to  50  km.  are  in  agreement.  While 
more  significant  differences  can  be  noted  at  ranges  closer  to  the  source,  if  should  be  noted  that  these 
differences  are  primarily  due  to  the  differences  in  how  PE  and  normal  mode  models  treat  a  layered 
environment.  Other  comparisons  (not  shown)  with  a  PE  based  spectral  decompo.sition  technique 
(Norton  and  Keiffer,  1991)  for  the  same  environment  show  very  good  agreement  at  all  ranges.  In 
general,  both  techniques  show  an  initial  rapid  fall  off  in  reflected  energy  then  a  leveling  off 
followed  by  an  increase  near  where  the  convergence  zone  occurs.  In  the  convergence  zone  the 
agreement  is  panicularly  good.  Clearly  the  main  features  in  these  predictions  correspond  quite 
closely  to  the  near-surface  propagation  indicated  in  Fig.  1. 


Fig.  2  Comparison  of  total  surface-reflected  intensity  from  PWS  extraction 
technique  and  from  normal  mode  theory. 


Now  that  it  has  been  verified  that  the  surface-reflected  PWS  extraction  technique  yields  the  correa 
reflected  intensity  as  a  function  of  range,  ray  theory  is  employed  to  verify  that  the  angular 
distribution  of  the  reflected  energy  as  a  function  of  range  is  also  correct.  There  are  three  regions 
that  can  be  identified  from  Fig.  1  that  can  be  clearly  associated  with  different  dominant  propagation 
paths.  Region  1  extends  from  the  source  to  approximately  5  km  and  is  identified  with  direct  paths 
to  the  surface.  In  this  region  the  PWS  extraction  techniques  predicts  essentially  equal  contributions 
to  the  reflected  intensity  from  all  angles  ve^  close  to  the  source,  by  2  km  the  contribution  is  limited 
to  15  deg.  and  bottom  bounce  paths  appeauing  at  very  high  angles  begin  to  become  significant.  At  a 
range  of  5  km  the  direct  path  is  limited  to  3  -  4  deg.  and  immediately  beyond  5  km  this  path 
abruptly  cuts  off  (precisely  as  predicted  by  ray  theory).  This  begins  Region  2  which  is  characterized 
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by  the  domuuuKe  of  bottom  bounce  paths.  As  the  rtinge  is  incmased  the  bosts','!!  boiaue  p.i-hs  unrie 
in  at  progressively  shallower  angles  (again  in  agreement  wnh  ihe  ra>  traces  i  untsl  about  40  pro 
where  the  convergence  zone  begins.  In  this  third  region,  two  paths  dominate  rays  initial!)  launched 
downward  at  shallow  angles  that  just  miss  the  bottom  and  those  that  undergo  some  inieractKin  wtt.h 
the  bottom  In  addition  there  are  rays  that  turn  very  near  the  surface  m  this  region  Ii  is  pre».ise;v  bic 
need  to  distinguish  this  latter  set  of  rays  from  those  that  actually  do  mteract  with  the  surf  ace  that 
makes  looking  at  PWS  at  two  range  separated  shon  anay  an  attractive  idea 

SUMMARY 

A  technique  has  been  described  that  can  extract  from  certain  underwater  propagation  models  an 
estimate  of  the  acoustic  field  incident  on  the  sea  surface  in  terms  of  the  PWS.  Thus  rs  useful  as  mpu! 
to  scattering  models  and  thus  comprises  one  pan  of  a  coupling  scheme  that  would  a!!nw  tor  the 
generation  of  hybrid  propagation/scattering  models  In  benchmarking  tests,  the  PWS  extraction 
technique  is  shown  to  provide  both  the  correct  to  reflected  energy  and  the  conect  angular 
dependence. 
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